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Three-wire thermocouple: Frequency response in constant flow 
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Theory and experimental measurements are compared with a novel three-wire thermocouple. 

Signals from three wires of unequal diameters are recorded from the thermocouple suspended in 
constant flow with a periodic temperature fluctuation. It is demonstrated that the reconstructed 
signal from the three-wire thermocouple requires no compensation for (x)^5(o lf where a > : is the 
natural frequency of the smaller wire. The latter result represents a significant improvement 
compared to previous work with two-wire thermocouples. A correction factor has also been derived 
to account for wires of arbitrary diameter. © 1995 American Institute of Physics. 


I. INTRODUCTION 

The evaluation of jet engine performance and fundamen- 
tal studies of combustion phenomena depend on the mea- 
surement of turbulent fluctuating temperatures of the gas 
within the engine. Historically, these temperatures have been 
measured with thermocouples. 1 The advantages of thermo- 
couples are their low cost, reliability, and simplicity since 
they do not require optical access or elaborate support elec- 
tronics. However, the design of a thermocouple represents a 
compromise between accuracy, ruggedness, and rapidity of 
response. 

For example, the measurement of fluctuating tempera- 
tures in the high-speed exhaust of a gas turbine engine com- 
bustor is required to characterize the local gas density gradi- 
ents or convective heat transfer. 2 Although thermocouples 
are suitable for the measurement of high-frequency tempera- 
ture fluctuations (<1 kHz) in a flowing gas or liquid, the 
measured signal must be compensated since the frequency of 
the time-dependent fluid temperature is normally much 
higher than the natural frequency of the thermocouple 
probe. 3 Moreover, use of a single-wire thermocouple in con- 
stant velocity flows requires knowledge of the fluid velocity 
and properties (e.g., viscosity, density, etc.) to determine the 
natural frequency. 

The present paper describes the performance of a novel 
three-wire thermocouple of unequal diameters that does not 
require compensation at lower fluid temperature frequencies 
nor any knowledge of the fluid velocity or properties. 4 5 The 
results of experimental measurements are presented along 
with the suggested procedure for the reduction of the data 
from the three-wire thermocouple as shown in Fig. 1. 


II. THEORY 

Use of a single-wire thermocouple requires knowledge 
of the fluid velocity and properties to determine the natural 
frequency. The latter quantity is necessary to establish a fre- 
quency dependent compensation factor for the measured sig- 
nal. 


A. Natural frequency of single wire 

The conservation of energy for a single wire is 
dT 

— ( 1 ) 


where the natural frequency 


4/t 

pcD * 


( 2 ) 


Here <o n is the natural frequency (s -1 ), T the temperature of 
the thermocouple (K), T g the temperature of the gas (K), t 
the time (s), p the wire density (kg/m 3 ), c the heat capacity of 
wire (J/kg K), h the heat transfer coefficient (W/m 2 K), and 
D the wire diameter (m). 

It is convenient to rewrite the heat transfer coefficient in 
terms of the Nusselt number 


h=^ Nu, (3) 

where k g is the thermal conductivity of the gas (W/m K) and 
Nu the Nusselt number. Thus the natural frequency in Eq. (2) 
becomes 


Ak g Nu 

where 6 

Nu=C 0 Re m Pr 1/3 

and 


(4) 

(5) 
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Here, Re and Pr are the Reynolds and Prandtl numbers, re- 
spectively, and ft is the viscosity of gas (kg/m s), v g the 
kinematic viscosity of gas (m 2 s _l ), and V the velocity of gas 
(m/s). 
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(8) 



FIG. 1. Schematic of three -wire thermocouple. 


The natural frequency of a single wire from Eqs. (4) and 
( 5 ) can now be written in the form 

k g C 0 4 Re m Pr 1/3 

72 . W 




( pc)D 2 

or separating the wire diameter 

a> n = CD m ~ 2 , 


(7) 


where the coefficient C depends on both fluid and wire pa- 
rameters, or 


4k e C 0 V m Pr m 

C= — — 

Substituting values corresponding to a type-K thermocouple 
and a gas velocity V=25 m/s for air at standard conditions 
into Eqs. (7) and ( 8 ), one obtains values of the natural fre- 
quency: 


<w„ (s *) 

D (m) 

D (mil) 

52.6 

50.8X10" 6 

2 

18.3 

101.6X10 -6 

4 

6.5 

203.2X10 -6 

8 


B. Determination of gas temperature 

The use of a single-wire thermocouple requires knowl- 
edge of the fluid velocity and additional properties such as 
viscosity or density. There are three unknowns that appear in 
Eqs. (1) and (7): (a) the coefficient C given by Eq. ( 8 ), (b) m, 
the exponent of the wire Reynolds number in Eqs. ( 6 )-( 8 ), 
and (c) the gas temperature T g . 

Consider three thermocouples of unequal diameters 
D l <D 2 <D 3 as shown in Fig. 1. The conservation of energy 
equation for each wire reduces to 


dT x , 

-^=CD”~ 2 {T g -T 1 ), 

(9) 

dT 2 

-^=CD?~ 2 (T g -T 2 ), 

(10) 

dT 3 7 


— = CD?- 2 (T g -T 3 ). 

(ID 


In principle it is possible to solve for the three unknowns 
that appear in Eqs. (9)— (11) in terms of the instantaneous 
values of the wire temperatures T j, T 2y and T 3 and their 
derivatives. However, because of the empirical nature of the 
Nusselt number given in Eq. (5) and uncertainty in the wire 
diameters near the thermocouple junctions, we choose to 
eliminate the coefficient C and solve for the gas temperature 
T g from Eqs. (9)— (11). Moreover, because the conservation 
of energy equations are linear at constant ambient velocities, 
we choose to use transformed data profiles. 

Thus taking the fast Fourier transform (FFT) of Eqs. 
(9)-(ll), squaring Eq. (10), and dividing by the product of 
Eqs. (9) and (11), one obtains 


( fl ) 

_/ D\ p 2 

(T g -T 2 ) 2 

\tjj 

\d x dJ 

(f g -r,)(r g -7',) 


Substituting 

r = t\lT- i t \ (13) 


into Eq. (12) and anticipating that the wire diameter ratio 
D 2 2 /D x D 2 ~0(\) or 


D 2 

D x Dy 


= 1 + e. 


(14) 
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the transform of the gas temperature becomes 


TABLE I. Wire diameter combinations ( m — 112 ). 


T = 
1 8 


K M 4u 0 c 0 

2 a 0 


(15) 


Here, b 0 = b + 2 et 2 + 0(e) 2 , -4a 0 c 0 = 4eT 2 2 (r- 1) 
+ 0(f) 2 , ho-4a 0 c 0 = h 2 +ey+0(e) 2 , and a 0 = r-l-e, 
where 


b = 2f 2 -r{T 3 + t x ) 


(16) 


and 


y=4T 2 [b + T 2 (r-l)]. 

We now seek a Taylor-series expansion of Eq. (15) for 
the transformed gas temperature about the point e=0. Thus 
to first order in e where 

f(e)= \jbl~4a 0 c 0 = y]b 2 + ey 
and 



1 

r-l-e’ 


one obtains 


f(e) = b+j£ e+0(e) 2 
and 


(17) 


1 e 

g( e )=~+ -^+°(e) , (18) 

where 

a = r~ 1. (19) 

Substituting Eqs. (17) and (18) into Eq. (15) and retaining 
terms 0(e), one obtains 



where a and b are defined by Eqs. (13), (19), and (16), 
respectively. 

We now define the transform of the gas temperature for 
the particular case where the wire diameter ratio 
D 2 ID x D^= \ in the form 

T g W— b ~ (21) 


or 


7,(0) = 


f 2 (f 2 f, + ^3-27^3) 


( 22 ) 


Substituting for the quantity b defined by Eq. (21) into Eq. 
(20), one obtains the final result 


where 


a = 



T X T 3 


(23) 


D^mil) 

Z> 2 (mil) 

£» 3 (mil) 

€ 

2 

4 

8 

0.0 

2 

3 

5 

0.17 

2 

3 

4 

-0.16 


and 7,(0) is defined by Eq. (22). The inverse transform of 
T g (e) becomes the reconstructed gas temperature or 

7 ’* = FFT " 1 [^(*)] (24) 

for assumed values of e. 

It is useful to consider the effects of wire diameter on the 
magnitude of the perturbation parameter 6. In general it is 
desirable to have small diameter wires so that their natural 
frequencies are shifted to larger values. Listed in Table I are 
three possible wire combinations that are commercially 
available along with the magnitude of e. The values of e are 
computed from Eq. (14) assuming that the wire diameters are 
exactly as listed by the manufacturer, which is very unlikely 
in practice. For example, a typical value of the exponent 
m = l/2 in Eq. (14) suggests that a 7% error in D 2 would 
yield a value of e^0.2 for the wire diameters listed in the 
first row of Table I. 


III. EXPERIMENT 
A. Apparatus 

In the present experiment, thermocouple sensors are ex- 
posed to a constant velocity airstream (<25 m/s) of varying 
temperature. In particular, the dynamic response of the ther- 
mocouple is measured for a periodic temperature profile of 
varying frequency. A rotating wheel configuration is used to 
deliver the test airstream to the proposed sensors. A similar 
experimental apparatus was described in detail by Elmore 
et al? and Forney et al. s 

A schematic of the rotating wheel apparatus used in the 
present experiment is shown in Fig. 2. As the wheel rotates, 
holes pass the two air supply tubes (3/4-in. i.d. copper) that 
allow slugs of hot (-55 °C) and cold (-30 °C) air to alter- 
nately enter a transition tube assembly mounted directly 
above the rotating wheel. In the transition tube the slugs of 
hot and cold air coalesce into a single airstream providing a 
periodic temperature profile covering a range of gas tempera- 
ture frequencies from roughly 1 to 60 Hz. A 3/4 -1/2 in. 
smooth copper adapter was inserted at the end of the coa- 
lescing tube (location of thermocouple) in Fig. 2 to increase 
the air velocity to 25 m/s. 

The analog temperature signal is digitized with a Data 
Translation DT-2801 A/D board mounted in an expansion 
slot of an IBM AT compatible computer as shown in Fig. 2. 
ASYST software was loaded onto the hard disk of the per- 
sonal computer and this provided a flexible system for data 
storage, manipulation, and display. The true temperature pro- 
file of the airstream is measured with a constant current an- 
emometer (TSI 1054- A) and sensor (1210-T1.5). 

The three type-K thermocouple wires were of diameter 
D 1 = 50.8, (2 mils), D 2 = 101.6, and Z) 3 = 203.2 fx m as 
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FIG. 2. Schematic of rotating wheel apparatus. 


shown in Fig. 1. The thermocouple junctions were fabricated 
by the Research Instrumentation Branch at the NASA Lewis 
Center. The wires were cut with a razor blade to produce a 
flat edge perpendicular to the wire axis. The wire segments 
were then mounted in a fixture such that the two faces of the 
junction were held together by springs. The junction was 
then formed by laser heating. The laser used was a KORAD 
model KWD Nd:YAG laser, operating at a wavelength of 
1060 nra. The power settings depend on the wire size and 
material, but for the thermocouples described in this report, a 
pulse duration of 4 ms was used, delivering a total energy of 
approximately 2 J. 

B. Procedure 

ASYST software was developed to acquire temperature 
data sequentially from the thermocouples and constant cur- 
rent anemometer. The data are digitized for four channels at 
a sampling rate of 512 Hz per channel for a total sample time 
of 0.5 s. After data acquisition, the FFT is taken for each 
channel. The largest peak (first harmonic) in the amplitude of 
the FFT is located for each channel which provides the fre- 
quency of the temperature profile. The ASYST software also 
records the amplitude ratio of the first harmonic for the 50.8 
jxm thermocouple-to-anemometer output. The amplitude ra- 
tio was recorded for each frequency setting of the rotating 
wheel and provided a value for the natural frequency (i) { of 
the small wire. 

The FFT of the signal from each thermocouple channel 
was substituted into Eq. (19) by the ASYST code along with 
an estimate of the parameter e. The inverse transform of f g 


by the ASYST code provided a reconstruction of the true gas 
temperature. 


IV. RESULTS AND DISCUSSION 

The temperature profiles for the ambient gas (large am- 
plitude) and the three type-K thermocouple wires are shown 
in Fig. 3. The three wire diameters used in the present study 
are D!=50.8, (2 mils), D 2 = 101.6, and D 3 = 203.2 /mi. 
The angular frequency of the gas temperature in Fig. 3 is 
cu=86 s 1 ('-T4 Hz) while the natural frequency of the small 



FIG. 3. Temperature profiles for the gas (large amplitude) and three type-K 
thermocouple wires. Three wires: D^SO.8 /im, &>j=42 s \ D 2 = 101.6 
>um, D 3 = 203.2 /tm. Gas temperature frequency: ou=86 s -1 . 
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FIG. 4. Gas temperature profile compared to reconstructed temperature de- 
rived from three-wire thermocouple, = 61 s _1 , o>=209 €=— 0.2. 

wire 0^=42 s” 1 . The natural frequency of the small wire 
corresponds to an ambient gas velocity of roughly 10 m/s. 

Figures 4, 5, and 6 illustrate the effects of increasing the 
magnitude of the wire diameter parameter 6 defined by Eq. 
(14). The small phase angle between the large ambient gas 
temperature profiles that appear in Figs. 4 and 5 and the 
smaller temperature profiles reconstructed from the thermo- 
couple probe allows one to distinguish between both profiles 
in Fig. 6. In the latter case, the temperature profile on the left 
is that of the ambient gas while the profile on the right was 
determined by taking the inverse FFT of T g (e) defined by 
Eq. (23) with 6=0.2. In all three figures the natural fre- 
quency of the small wire 6^=61 s -1 (V— 20 ms -1 ) while the 
angular frequency of the ambient gas temperature is <o=209 
s _1 (-33 Hz). 

It is clear that a value of the parameter 6=0.2 provides a 
signal that requires no compensation for the frequency ratio 
a)/ w 1 =3 .4. It should also be noted that a 7% error in the 
diameter of D 2 would account for a value of 6=0.2. Further, 
the coefficient of 6 in Eq. (23) depends on the difference 
between the uncompensated FFT of the signal from the wire 
of diameter D 2 and the compensated FFT of the signal from 
the three-wire probe with a particular diameter ratio 



TIME (SEC) 


FIG. 5. Same conditions as Fig. 4. e=0.0. 
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FIG. 6. Same conditions as Fig. 4. e=0.2. 

Dl=DiD 3 (e-0). Thus the contribution to the signal by the 
second term on the right of Eq, (23) is negligible at low 
frequencies but increases with larger values of oi/oij . 

The amplitude ratio of the three-wire thermocouple 
output-to-true gas temperature as a function of gas tempera- 
ture frequency has been plotted in Fig. 7. Also shown is the 
first-order response of the small wire. The open symbols cor- 
respond to a value of 6=0 in Eq. (23) where the amplitude 
ratio requires no compensation for a)^2w l . For the range of 
values 2w 1 <o><5w 1 the amplitude ratio follows the function 

- 4 = i - 92 Uj ■ < 2 « 

It should be noted that the first-order response is also propor- 
tional to (cu/ct)!) -1 at large frequencies. 

A significant improvement is observed, however, when 
the parameter 6=0.2 as indicated by the solid symbols in 
Fig. 7. For the latter case, no compensation is required for 



co/co 1 


FIG. 7. Amplitude ratio of three-wire thermocouple output-to-true gas tem- 
perature. Aiso shown is first-order response of small wire. 
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(0^5 (Di . As discussed earlier, the amplitude ratio data co- 
elesce at low gas temperature frequencies where the coeffi- 
cient of € in Eq. (23) approaches zero. 

Both values of e=0.0 or 0.2 are of little value, however, 
when (o^6a) x where the amplitude ratio for the three- wire 
thermocouple reduces to the first-order response of the 
smaller wire. The latter result appears to occur at a gas fre- 
quency of —60 Hz. For the present experimental system the 
signal-to-noise ratio may be responsible for the small re- 
sponse at large frequencies. This situation may be improved 
by replacing the two larger wires of 4 and 8 mils diameter by 
wires of 3 and 4 mils, respectively, as listed in Table I. 
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